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I. INTRODUCTION 

By the recent observations of the universe, we have 
found that the universe is undergoing a phase of acceler- 
ated expansion at the present epoch from about 5 x 10^*^ 
years ago and the universe is flat. 

The accelerated expansion of the universe has been es- 
tablished by the observation of the type la Supernovae 
(SNIa) HI. This type of supernovae has an abosrption 
spectrum of Silicon, which helps us to distinguish the su- 
pernovae from other types. An important thing is that 
the intrinsic luminosities of SNIa are almost uniform and 
their absolute maginitude is almost —20, from which we 
can find the distance between the earth and the super- 
novae. From the distance, we can find how old the SNe 
are. We can also the speed of the SNe from the redshift. 
Combining the distance and the redshift, the accelerated 
expansion of the universe has been established. 

The flatness of the universe can be found from 
the observation of the anisotropics of Cosmic Mi- 
crowave Background (CMB) observed by the baloons, the 
BOOMERANG project 3 and the MAXIMA-1 project 
. The spectrum of the CMB has a characteristic struc- 
ture and the CMB was produced in the early universe, 
that is, when the universe becomes transparent or the 
mean free path of the light becomes infinite. Then the 
CMB has propagated very long distance, about 10 bil- 
lion light years and if the univers is curved, the structure 
of the spectrum should be deformed. The deformation 
has, however, not observed, which tells that the universe 
should be flat. 

If the universe (the spacial part) is flat, we may assume 
the metric of the flat FRW universe in the following form: 

ds'' ^ -dt^ + a{tf {dx'f (1) 

1=1,2,3 

Here a is called a scale factor. If a > 0, a > 0, the 
universe is expanding and accelerating. Related with the 
accelerated expansion, we often use a parameter w, which 
is called as "equation of state" parameter. The present 



data tell 

w=-^-l. (2) 
P 

Here p is the pressure and p is the energy density of the 
universe. We may explain the relation between w and the 
accelerated expansion later but if u; < — I, the universe 
is accelerating, which Eq.(|2Jl tells. As w is negative, if 
the energy density p is positive as usual, the pressure p 
should be negative. 

In case of the radiation, we have w = ^ since the trace 
of energy-momentum tensor vanishes: T = —p -I- 3p = 0. 
Usual matters, like baryon, can be regarded as dust with 
w = (p = 0). The cosmological term can gives w = —\. 
It is, however unnatural if the cosmological term, which 
may be suggested from (jJJ, is a part of the fundamental 
gravity. The cosmological constant should be given by 
A ^ i?^, where H is the Hubble parameter, H = ^. Since 
the Hubble parameter is given by ^ 70 kms^^Mpc"^ 
~ 10^^'^ eV, the scale of the cosmological constant is 
unrealistically very small compared with the fundamental 
scale of the gravity, the Planck scale 10^^ GeV = 10^* eV. 

On the other hand, since the universe is flat, the energy 
density of the universe should be almost critical density : 
10~2^ g/cm'^. The density of the non-relativistic matter 
(baryons, dark matter) is, however, almost | of the crit- 
ical density. This tells that 70 % of the energy density 
of the universe should come from the matter (?) with 
negative w. We call the matter (?) as dark energy. Con- 
versely, the density of the dark energy is almost same 
order with the density of the usual matter. Then maybe 
it would be more natural to consider that the dark energy 
would be a kind of matter, rather than that it would come 
from the fundamental cosmological term in the gravity. 

Now we explain the relation between the accelerated 
expansion and w. The energy conservation law = 
V^Tfj^i, has the following form in the FRW metric |^: 

p + ?,H{p + p)^ p + i{l + w)Hp . (3) 

In case that is a constant, since H = we can inte- 
grate (|3j) : 

p = poa-3(i+-). 
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The (1st) FRW equation, which is the (t,t) component 



2 



of the Einstein equation, is given by 
6 , 



(5) 



Here = IGttG and G is the gravitational constant. 
Then if w ^ —1, the scale factor a can be found by 
integrating (0) by using 



(6) 



The case w ^ —1 corresponds to deSitter universe: 



floe 



(7) 



Then if it; < 



that is, 



> 1, the universe is 



3{w+l) 

accelerating, which can be found more directly, by com- 
bining 2nd FRW equation, 

"2 2 

- = -^ip + 3p) = -^{l + 3w)p , (8) 
a 12 12 



which tells d > if w < 

If— 1 < w < — ^, the universe is expanding and 
accelerating and the corresponding matter is called as 
"quintessence" 0|. If w < — 1, the universe is accelerat- 
ing but shrinking in © but if we change the direction of 
time as t —^ ts — t, the universe is expanding (and acceler- 
ating). The corresponding matter is called as "phantom" 
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As clear from after replacing t ^ tg — t, when 
w < —1, there is a singularity at t — ts, which is called 
the Big Rip singularity. We will mention more about the 
singularity later. The present data of the universe seems 
to tell that w might be less than —1. In the viewpoint of 
the quantum field theory, however, this is rather strange. 
As an example, we consider the scalar field with the 
potential V. If we assume the FRW metric and (p only 
depends on the time coordinate, the enrgy density p and 
the pressure p are given by 



1 • 



1 



(9) 



Then if V{(f>) > ^<p^ , we have negative w but the value 
of w should be greater than —1. In order that u; < — 1, 
we should change the sign of the kinetic energy: 

P^~\^' + V{^), p^A^^^V{cl>), (10) 

which tells that the phantom might be ghost with neg- 
ative norm, which should be a serious problem in the 
quantum field theory |3ll |. Conversely, if surely w < —1, 
there should be a new physics. 



II. SCALAR WITH EXPONENTIAL 
POTENTIAL 

In this section, as a most simple model, we consider the 
scalar field theory with the exponential potential. This 



model has been investigated from a very long time ago, 
say in 't'I . The action of the model is given by 



S = 



F(0) 



Vbe *o 



(11) 



Usually 7 > 0, but for phantom, we may have 7 < 0. We 
assume (flat) FRW universe and (p only depend on 
time coordinate' t. Then the ^-equation fo motion and 
the FRW equation has the following form: 







-7 



(12) 



If we assume </) oc In t, /i cx i -'^ , we can solve the equations 
in H12II very easily : 



.In 



H 



At 



Then the scalae factor a is given by 

t 

a = ao — 



Since a = aot2(™+'^' for general w, we find 

-1 



w 



(13) 



(14) 



(15) 



Then if 7 < 0, surely phantom appears. 

A very interesting point is that the general solutions of 
the scalar model with exponential potential can be found. 
In case 7 > case 0, if we assume 



i±i , 2(w - u) 
a = e 3 , (b = ■ 



3Vq 

dr — dt\l e <t>o^ , 



(16) 



the Hamiltonian constraint and other equations have the 
following form: 



and 



Here 



dv du ^ 
dr dr 



I, -2^„ c''"^ I, -2\„ 



u 



(17) 



(18) 



(19) 
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The equations in H17() can be solved very easily and when in the Einstein frame one: 



a < 1, the behavior at t — > +00 coincides with that of 
the previous solution H13|l . H14|l and w has the following 
form: 



(20) 



Then if < |, the universe is expanding and accelerat- 
ing (quintessence). 

The case 7 < has been also exactly solved in 0, 
instead of H16|l. by using a complex variable z, we assume 



e 3 



K0o\/-37 



Z = e^i-"^)^ 



(21) 



Then the Hamiltonian constraint and other equations 
have the following form: 



dz dz* ^ 
dr dr ' 



d^Z 
1^ 



= (1 



(22) 



The above equations can be also solved easily. When 
i ^ (if we replace t ^ ts — t, t ^ tg), the behaviour 
corresponds to the one in the previous solution ad we find 
w is given by 



w 



1 -2 1 

-I - a < -1 , 



(23) 



which surely corresponds to the phantom. 



III. A MODIFICATION OF EXPONENTIAL 
POTENTIAL MODEL 

Since the models with negative norm would not be con- 
sistent with quantum field theory as the theory may be 
non-unitary. We consider a modification of the model 
based on (9|. 

We may start with Jordan (string) frame action, which 
is the Brans-Dicke type 10], 



S ^ ^ j d'^xy^c°"t' (i? 



+ I d''xV^{--d,xd''x-Uix) 



(24) 



Here a and 7 are constant parameters and 7 can be neg- 
ative. The scalar field x represents the matter. Since x 
does not interact directly with </>, the equivalence princi- 
ple would not be violated although the effective gravita- 
tional constant Ke^~ depends on (f>. 

By using the scale transformation : Qf^^ = 
e~'^^'^gE ^1}^, the Jordan frame action H24|) is transformed 



e-^'^y(0) 



7 A ^t.ug 



2 J 9e a^^cpOi, 

+ j d'^xy/^ 

-e-^^U{x) 
Even if 7 < 0, when 

{d-l)a^ 
2 



9e df^xduX 



> 



(25) 



(26) 



the kinetic energy becomes positive. Especially when d — 
4 and 7 = —1, we have > ^. Then the sign of the 
kinetic energy depends on the frame. 

When d = 4, X = and V{(j)) is exponential type, we 
can solve the system in a way similar to the last section. 
If we define, 



if : 
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y(0) = Vbe ™ (27) 



and assume the FRW metric in Einstein frame : 



-dtl 



(28) 



1=1,2,3 



we obtain the following solution 



3 ,„2 
4V0 



ttE = aEo [ - — , ip^ ipoin - — 

^tEO / tEO 



1 /27 , 3 



(29) 



For the FRW metric in the original Jordan frame in 
we have 



dt = e'-t'dtE 



(30) 



Then we obtain 



a = e ^ ttE = cleo 



t \ i-iM) 



a 1 



Inl 

Pvo \ to 



(■ »y 

Then the condition for the acceleration is given by 



(31) 



1 - 



>1 , 



(32) 
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and w is given by 



3,^2 Vpa 



(PO {2(3 - 9v?o) + 8 
3 (2/3 - 3(^o) "^0 



(33) 



(34) 



Since the denominator in the above expression can van- 
ish, by properly choosing the parameters a, ipo, (7), we 
can put w in an arbitrary value. Especially if = 

2, we have w = —1 and when (^q > | we find 

^ < -1 if 1^0 > ^f§= > 2 
->-l if -^>|;§or-^<2. 

On the other hand, when yjg < |, 
Even if 7 > 0, we can have w < —1. For example, if 



(35) 



^0 = 4, ^=3(/3.-L 

we have w, which is less than —1: 
11a/2 + 203 



(36) 



213 



= -1.025... < -1 . (37) 



IV. l/i?-MODEL AND GENERALIZATION 

In a modification of Einstein gravity has been 

proposed. In this section, we consider the model and its 
generalization (see also [13. Il^ l. 

The action of the model in (lH include the inverse 
power of the curvature: 



By the variation of B, we obtain A = R. By substituting 
this equation into the above action (|4U|I . we obtain the 
original action. If we consider the variation of A first, 
we obtain B = f'{A), which can be solved as A ~ g{B). 
Then we obtain 



S 



1 



d^x^g{B{R-g{B)) + f{g{B))}. (41) 



Instead of A, we may delete B as follows, 
1 



5 = 



d'xV^g{f{A){R-A) + f{A)}, (42) 



which can be regarded as the Jordan frame action. By 
using the scale transformation g^;^ giiv with a = 

— In /'(A), we obtain the Einstein frame action : 



Se — 



1 



(Tx^/^i R 



9 d.Ad.A f'{A) + nAr 



f'iA) 



1 

V{a) = e' 
A 



^[R- "-gP^dpad^a - V{a) ) (43) 

2(j _ 



d^X 

5(e-'^)-e-/(g(e-)) 



nA) nAf 



(44) 



Especially, in case of ll'l, when A(= K) is small, V[a) 
behaves as a exponential function: 



V{<j) 



re3^ 



(45) 



Then we can solve the equations as in the previous sec- 
tion. Even for more general case [15| as 



/(A)=^ + 7A-"(ln^ 



(46) 



when curvature is small and n 7^ 0, the potential behaves 
as an exponential function 



S- 



d^Xy/^ R - 



R 



(38) 



When curvature is small {\R\ <C /i^), we have a (x and 
w = — |. This model may be called as "c-essence" ("c" 
expresses the curvature) [l4j . 

We may consider the generalization of the model as 
follows, 



(39) 



Here /(i?) can be an arbitrary function. By introducing 
the auxilliary fields A, B, we may rewrite the action 
in the following form: 



d^xV^{B{R-A) + f{A)} 



(40) 



.n+1 

and we can solve again the equations: 



a ^ t~ 



w 



3(n-f l)(2n+ 1) 



Then we find that when 

-7-6V2 



n < 

or n > 
w < and when 



12 

-7 + 6\/2 



l<n<-- 



12 



n < 



l + \/3 



or n > 



-1 + V3 



(47) 



(48) 



(49) 



(50) 
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the universe accelerates. 

In case that there is no matter and the Ricci tensor 
R^,y is covariantly constant (Vpi?^^ = 0, which means 
Rfiv = const. X g^,^), the equation of motion is reduced 
to 



= 2/(i?) - Rf{R) , 



(51) 



which is an algebraic equation with respect to R. For 
example, when 



fiR) = R 



we have 



-R 



(n + 2)a 



+ (m - 2)bR" 



(52) 



(53) 



Especially, in case of {n — 1, a — fi"^ , b — 0), we find 
R = ±V3^M^ . (54) 

The + sign corresponds to deSitter space, which is an 
exponentially expanding and accelerating universe. We 
may regard this solution would correspond to the present 
universe. We should not that even if 6 7^ but m ~ 2, 
we obtain the same solution as (|54|l . 

Just after the l/i?-model was proposed, there were sev- 
eral criticism [T^ and an improvement was proposed 
in ^3 J where the action is modified as 



f(R)^R- +bR^ 
It 



(55) 



Simultaneously with the acceleration of the present uni- 
verse, the model may explain the inflation at the early 
stage. The models of the inflation using i?^-theory has 
already been investigated in 0, We should note 

that when R is large, where n goes to 2, from H47|l . we 
find w ^ — 1, which corresponds to the deSitter space. 

First, we consider the instability pointed out in [T6l| . 
The general equation of motion with matter in f{R)- 
gravity is given by 

\9^.uf{.R) - R^.uf{R) - g^unfiR) + V^V./'(i?) 



-T, 



(56) 



Here T(^m)fii' is the matter energy-momentum tensor. By 
multipling the above equation with g^'^ , we obtain 



DR 



,f(^){R) " ^ 3/(2)(i?) 3/(2)(i?) 



6/(2) (i?) 



T 



(57) 



Here T = ■, We consider a perturbation from the 
solution of the Einstein gravity: 

^2 

R = Ro = -—T > . 



We should note that T is negative since \p\ ^ p in the 
earth and T = —p-\-Sp^ —p. Then we assume 



R 

Then we find 







DRo 



i?o + i?i , (|i?i| « |i?o|) . (59) 

f'{Ro)Ro 



2/(i?o) 



VpRoV^Ro 



Rn 



3/(2) (i?o) 



3/(2) (i?o) 3/(2) (i?o) 

j^'Hro) 



-DRi 



VpRoVRi + UiRo)Ri 



UiRo) 



/(2)(i?o 
^v/(2)(i?o) f^'HRor 

f('HR^)f^^\Ro)Ro fW(Ro) 
3/(2)(i?o)2 3/(2) (i?o) 

, 2f{Ro)f ^\Ro) i?o/(^) 



VpRo^PRo + -Ro 



3/(2)(i?o)2 3/(2)(i?o)2 



(60) 



If U{Ro) is positive, since DRi ^ —dfRi, the pertur- 
bation Ri becomes exponentially large and the system 
becomes unstable. We may regard Vpi?o ~ if we as- 
sume the matter is almost uniform as inside the earth. 
In case of the model in we find 



U{Ro) ^-Ro + f^ 



(61) 



If fi corresponds to the acceleration of the present uni- 
verse as in (|54f) . we obtain 



lO^^sec ' 



(lO-33eV)"' 



(62) 



This corresponds to the deSitter space which expands 
exponentialy. If the present universe expands with power 
law as in H47(l . the value of a can be much larger than 
the above value. 

The value of Rq has been evaluated in [T^ : 



^ ~ (l0-26s 

a ^ 



Pn 



Rq 



'lO^sec) 



gem 
P 



-3 



gem 



-3 



(63) 



26 



Then the system is unstable and would decay in 10 
sec. 

In the model [l3| with a = p^, if b :$> t-S^, we 
find 



U{Ro) - ^ > . 



(64) 



The system is unstable again but the decay time becomes 
about 1000 second and macroscopic, that is, it is im- 
proved by 102^. But the assumption b 3> is not so re- 

(58) alistic since this means b-^ < (lO"eV)^ = (l02GeV)^ 
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In "2^, it has been given the constraint for inflation in 
i?2_theory from COBE-DMR data of CMBR (cosmic mi- 
crowave background radiation) as 



M 
rripi 



2.6 X 10"^ 



(65) 



that is ~ (lO^^GeV) , which contradicts with the 
assumption. Anyway instabiUty depends on the details 
of f{R). Furthermore, if the present universe expands by 
the power law, we need not to assume 



lO^^sec . 



(lO^^^eV)" 



(66) 



In it has been pointed out the problem in the 
equivalence principle. By the scale transformation, the 
matter action is transformed as 



(67) 



Then the matter couples with cr, which may violate the 
equivalence principle due to the force mediated by a. In 
case of the model in if we assume H66|> again, the 
mass of (T becomes 10"^'^ eV= 10"''^ GeV and very light, 
which will violate the equivalence principle. 

In the model in the neighbourhood oi R = A = 
■\/3a — V^ii^ , we obtain 



d^V{a) 



dl) dA^ J 



(-i + 6\/3^) (| + 26V^)^ 



Then if 



1 



3v 3a 



A=V3a 



(68) 



(69) 



which means f'{A) ^ 0, the mass of cr becomes very large 
and would not conflict with the equivalence principle. In 
general, since 



Via) = e-5 (e-^) - e'^ f 



[9 e 



1) 



A 



r{A) /'{Ay 



(70) 



the mass of a would become large when f'{A) ^ 0. Then 
if the present acceleration of universe corresponds to de- 
Sitter solution i? = > 0, if we consider a model where 
/'{Aq) ~ 0, the mass of a can be large. 



V. BIG RIP SINGULARITY 

In case of phantom (w < —1), the scale factor of the 

2 

universe behaves as a = ag (i^ — t) . Then a di- 
verges at t ^ ts- Then anything like atom or hadron. 



which has internal structure, should be torn. Then this 
singularity is called the Big Rip singularity ( "rip" meas a 
long tear or cut). This is the third possibility of the fate 
in the universe, following the Big Crunch and the eternal 
exapnsion. 

We now consider the origin of the Big Rip. First we 
note that the energy density p behaves as 

-3(H-u;) 



P = Poa 



(71) 



Then in case of phantom (w < —1), when a increases, p 
also increases. 

In order to consider any concrete model, we consider 
the scalar field theory coupled with gravity, whose action 
is given by 

S^-^l d^x^g (r - l^d^cfd^cj, ~ V{4>)) . (72) 

If 7 < 0, we have a phantom. By combining the equation 
given by the variation of 0: 



= -7 ( ^ + 3iJ 



y'(0) 



and the (1st) FRW equation : 



1 f dip 



we find 



(73) 



(74) 



(75) 



which is positive if 7 < 0, > 0, and 7^ 0. Then the 
energy density of phantom (7 < 0) increases with time. 

We now consider if the Big Rip singularity really exist 
or not. First we should note that since 7 < 0, we have 
P41 < ^(0)- Then if there is an upper bound in V{(j)), 
there is also an upper bound in p^. The potential V{(l>) 
might behave exponentially in the present epoch of the 
universe, but if there is an upper bound in V{(f>), there 
does not occur the Big Rip singularity and the space 
becomes deSitter space asymptotically. 

Another possibility to avoid the Big Ri p s ingularity is 
to include the quantum correction [21| [32|- Near the 
Big Rip singularity, since a blows up, cuvatures become 
large as i? cx |i — <s| ^- Since the quantum correction 
contains power of the curvatures in general, the quantum 
correction becomes very important. In order to consider 
such a quantum correction, we include the contribution 
from conformal anomaly : 



T = b[ F+-DR 



b'G + h"UR , 



(76) 



Here F is the square of 4d Weyl tensor and G is the 
Gauss-Bonnet invariant: 



F — ^ 2i?y i?'-' -I- RijkiR^^ 

kl 



kl 



a 



R^ - AR,jR'^ + R,jkiR'' 



(77) 
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In case that there are N scalar, N1/2 Dirac spinor, Ni 
vector fields, N2 (= or 1) gravitons, ^Yhd higher deriva- 
tive confomial scalars, the coefficients b, and b" are 
given by 



TV + 6N1/2 + 127Vi + 6117V2 - 8A^HD 
120(47r)2 ' 
N + llA^i/2 + 62iVi + 1411iV2 - 287Vhd 



b' 
b" 



360(47r)2 



(78) 



We should note & > and 6' < for the usual matter ex- 
cept the higher derivative conformal scalars. We should 
note that b" can be shifted by the finite renormalization 
of the local counterterm . Then b" can be arbitrary. 
By using the corresponding energy density pA and pres- 
sure p^, Ta is given by Ta = — + 3pA- Then by using 
the energy conservation in the FRW universe 



= ^+3i7(pA+PA) 
at 



we may delete ^a as 



(79) 



Then we also have an expression for 'pA 
1 dp A 



PA 



'PA 



^ - 4 / dta-'HTA 



Ta 
3 

1 
3 



3H dt 
1 



, 'dH\ 



/2 „\ {d^H d^H 



dt^ 



\ I dta^H 
a* 



12i/^ 



,dH 
~dt 



-126 



-K246' ■ 



dH\ 



dt I 



dH 
'dt 

dH 
'dt 



-6(^6 + 6" 



d^H d^H 



[ dt^ 



,dH 
~dt 



(82) 



By including the quantum correction in (|81|l . we may 
write the corrected FRW equation in the following form: 



6 



7 



(83) 



Ta = -ApA - 



1 dpA 
H dt 



(80) 



We now assume 

H = ho + Sh , = 00 In 



t,-t 



(84) 



which gives the following expression of pA- 



PA 



X I dta^HTA 
a* 

' dta^H 



-126 



dH 
~dt 



\ dt^ 



dt^ 



dH 
~dt 



12H^ 



,dH 
'dt 



and when t tg, 6ti, S(f> are much smaller than the first 
terms but might be singular. Then (^-equation of 
motion reduces as 

= -7 



2Votl 



3/io 

t.~t 



^0 {ts ~ ty 

which gives 

Votj 



1 - 



+o((i,-t) 1) ,(85) 



1^1 



[ts - t) 



If we assume |6 -I- 6" 7^ 0, we find 



PA-^h[-b + b"]-^ 



(86) 



(87) 
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Then by substituting H86|) and H87() into the corrected 
FRW equation (|S3J), wc find 







and 



dh 



2 (§6 + 6") 



(ts - t) hi 



t2 



(88) 



(89) 



Here t2 is a constant of the integration. Then the scale 
factor a behaves as 



a — 



t,-t 



t2 



l(§b + 6"j 



xe . (90) 



Big Rip singularity, the quantum correction would be 
very important. 

There are many other models, as tachyon (see |23| 
for an example), the Born-Infeld theory (see_[2^ for 
example), Chameleon model l2al (see also [23 ), and 
R + imattcr-R" type modcl in [23. There are applica- 
tions of the modified gravities to the cosmology (for ex- 
ample, ll^). These theories have structures as in the 
models obtained by the dimensional reduction from the 
model in the higher dimension. In fact, several models 
are proposed from the viewpoint of the brane world (for 
example, see jsj)- Near w = —1, we feel that theories 
in the extra or higher dimensions (brane world, Kaluza- 
Klein model) would appear. 



There appear logarithmic singularities in ^ but 
the singularity becomes rather mild. Then the universe 
might develope beyond t — tg. 
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